End  (1) 


AD-A197  498 


ftWi  Dmt  £nitttbQ _ 

ITATiriU  READ  INSTRUCTIONS 

1 1  A  1 IUW  r  AO C _ BEFORE  COMPLETING  FORK 

I*.  OOVT  ACCESSION  NO.  L  RECIPIENT'S  CATALOG  HUMBER 


4.  TITLE  fan d  Submit) 

Advances  in  che  p  and  h-p  Versions  of  the 
Finite  Element  Method.  A  Survey 


S.  TYPE  OE  REPORT  *  PERIOD  COVEREO 


iii,: 


«-  performing  oro.  report  number 


7.  AUTMORf#) 


I.  Babuska 


«.  CONTRACT  OR  GRANT  HUMBERT  •) 

ONR  NOOO14-850K4169 
NSF  DMS-84-10324 


».  performing  organization  name  ano  address 

Institute  for  Physical  Science  and  Technology 

University  of  Maryland 

College  Park.  MD  20742  _  _ 


ii.  controlling  office  name  ano  aooress 

Department  of  the  Navy 
Office  of  Naval  Research  I  '*■  humoer  of  paces 

Arlington.  VA  22217 _ I  12 _ ; _ 


4.  MONITORING  AGENCY  name  A  AOORESV"  dlt/erent  frM  Controlling  Of  tic*)  I  IS.  SECURITY  CLASS.  (et  till*  report) 


m 


I*.  DISTRIBUTION  STATEMENT  (ol  ihlt  Rtperl) 


Approved  for  public  release:  distribution  unlimited 


*7.  DISTRIBUTION  STATEMENT  (et  the  eh*  tract  entered  M  Block  70,  II  different  item  Report) 


IS.  KEY  WOROS  (Continue  on  aide  If  nacaaaary  end  Identify  ky  block  number) 


20.  ABSTRACT  (Continue  on  reeorae  aide  It  neceaeery  end  Identity  by  block  number) 

The  paper  gives  the  survey  of  the  advances  in  the  theory  and  practice  of 
the  p  and  h-p  versions  of  the  finite  element  method.  It  gives  the 
extensive  list  of  references  related  to  recent  results  of  this  new 
approach. 


oo 


FORM 
JAN  71 


EDITION  OF  I  NOV  «S  IS  OBSOLETE 
5  N  0107-  IF-  014-  6401 


SECURITY  CLASSIFICATION  OF  THIS  PACE  (*>>tn  Obit  Snltrtb) 


} 


ADVANCES  IN  THE  p  AND  h-p  VERSIONS  OF  THE  FINITE  ELEMENT  METHOD.  A  SURVEY 


I .  Babuska 

Institute  for  Physical  Science  and  Technology 
University  of  Maryland,  College  Park,  MD  20742 

l)  n 

ABSTRACT .  The  paper  gives  the  survey  of  the  advances  in  the  theory  and 
practice  of  the  p  and  h-p  versions  of  the  finite  element  method.  It  gives 
the  extensive  list  of  references  related  to  recent  results  of  this  new 
approach . 

1 .  /  Introduction 

'“The  finite  element  rpethod  has  become  the  main  tool  in  computational  me¬ 
chanics.  The  MAKABASE</f42-j — f~4^t  contains  at  present  the  information  on  ap¬ 
proximately  1400  finite/boundary  element  programs,  about  20,000  references  on 
finite  element  and  2,000  boundary  element  technology.  To  date  there  are  more 
than  two  hundred  monographs  and  proceedings 7[44 ]  .<>  Recently  a  new  direction 
in  the  finite  element  theory  and  practice  appeared,  the  p  and  h-p  versions 
which  utilizes  high  degree  elements.  About  3-4  dozens  references  (out  of 
22,000)  and  only  few  programs  are  available.  The  aim  of  the  paper  is  to 
briefly  survey  the  state  of  the  art  about  the  p  and  h-p  versions  and  present 
the  basic  references^ 

2.  THE  MODEL  PROBLEM  AND  ITS  PROPERTIES  (•  .  . _ _ 

~  We  restrict  ourself  to  the  most  simple  b yt  e  characteristic  model  problem 
for  elliptic  partial  differential  equations.  Let  0  €  R  ,  ((x^,X2>  >  x)  be  a 

M 

simply  connected,  bounded  domain  with  the  boundary  30  -  T  -  U  T.;  are 

analytic  simple  arcs  called  edges  i-1 

-  {(<Pi(0,^1(0)U  €  I  -  [-1,1]} 

—  9  9 

where  (!)  ,  ^(!)  are  analytic  functions  on  I  and  |y>j(!)l  +  l^[(!)r 

>  >  0.  Let  A^,  i  -  1,...,M  be  vertices  of  O  and  -  A^A^+^  i.e. 

the  edge  is  linking  the  vertices  A^  and  A^+^.  By  i  -  1,...,M  we 

denote  the  internal  angles  of  O  at  A^. 

Let  us  be  interested  in  the  model  problem 


(2. 

1) 

-Au  -  f  on  01 

(2. 

2) 

u  -  >s>  on  r°  -  u  r.  =  0 

j€Q  J 
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Here  Q  is  a  subset  of  -  M.  The  performance  of  the  finite  ele¬ 

ment  method  strongly  depends  on  the  properties  of  the  (exact)  solution  of  the 
solved  problem  especially  on  its  smootheness . 

In  structural  mechanics  the  problems  are  typically  characterized  by 
piecewise  analytic  input  data.  In  our  case  the  boundary  is  piecewise  analytic 
and  we  will  also  assume  that  f  is  analytic  on  fi  and  <p  and  ^  are  ana¬ 
lytic  on  rr  The  question  arises,  how  to  best  characterize  the  solution  of 
the  problem  (2.1)  (2.2)  for  such  input  data.  Is  it  possible  to  do  in  the 
terms  of  countably  normed  spaces. 

Let  A1  -  i,x21),  r£  -  (xj/X^)2  +  (x2-x2)t)2.  P  -  {0^},  0  <  Pt 

JL  Pi 

<  1,  i  -  1,...,M.  Define  $a(x)  -  I  i  r  (x)  and  for  any  integer  k  let 

P  i-1  1 

M  0  +k  x 

90±kW  “  1  1  ri  (x).  Let  finally  8£  (A)  -  (u  €  HL(Q)  |  ||^+k  2Dau||  < 

Cd^k!  k  *  2,...,|a|  -  k,  C  and  d  independent  of  k).  We  denoted  Dau  - 

a\+a2  a2  o 

d  vi/dx  dx  ,  al  +  a2  "  I®!*  Function  u  €  B^(0)  is  analytic  on 
M  12 

fl  -  U  A,  and  has  a  specific  behavior  in  the  neighborhood  of  the  vertices 
i-1 

A^ ,  i  -  1 , . . . ,M.  Essential  is  that  the  solution  of  the  problems  of  struc¬ 
tural  mechanics  belongs  to  Bjg(fl) .  In  [8]  [11]  we  have  proven 

Theorem  2.1.  Let  f  be  analytic  on  0,  ip*  be  analytic  on  T  ,  j  € 

—0  —  J 

Q  and  continuous  on  T  ,  i/>  be  analytic  on  T  ,  j  €_M  -  Q.  Let  u  be  the 

solution  of  (2.1)  (2.2).  Then  u  €  8^(0)  witlr  0^  >  0^  where  depends 

on  and  the  type  of  boundary  conditions  on  D 

We  can  now  ask  the  question  how  to  completely  characterize  the  class  of 

the  functions  <p  and  ^  and  f  which  guarantees  that  the  solution  of  (2.1) 
(2.2)  belongs  to  8^(0).  This  has  been  done  in  [9],  Let  us  mention  this 
characterization  in  a  special  (simplified)  setting. 

O 

Theorem  2,2.  Let  Q  be  a  polygon,  u  e  8^(0)  and  g^  be  the  trace  of 
u  on  Then 

a)  if  1/2  <  0it0i+1  <  1.  then  gi  €  B*  (Tt) ,  €  (01+j_r  1/2, 1/2), 

j  “  » 2 • 

b)  if  0  <  PL,Pi+l  <  1/2,  then  gi  $  Bl  (i^),  €  (0i+j  _i+  1/2 , 1) . 

Here  we  denoted  analogously  as  before  for  s  -  1,2 


2  ^ 

8s <r.)  -  {g  €  Hs-1(ri) | UTT  Ix-a1+1.,  i^.:i+i'sg[*l 
P  j-i 


^L2(r.)  - 


In  [9]  we  also  have  shown  the  extension  of  the  traces  g  belonging  to 
the  spaces  mentioned  in  Theorem  2.2  to  the  function  from  8|(0).  By  this 


( 


characterization  and  Theorem  2.1  it  is  easy  to  verify  that  the  solution  of 
(2.1)  (2.2)  belongs  to  8£(0) . 

In  [14]  we  have  proven  that  the  eigenfunctions  belong  to  the  space 
81(0)  too. 

Although  we  have  discussed  only  the  special  problem  (2.1)  (2.2),  Theorem 
2.1  is  valid  under  much  more  general  assumption  and  also  for  system  of  equa¬ 
tions  [10].  We  also  have  characterized  the  countably  normed  spaces  which 
characterize  the  solution  of  (2.1)  (2.2)  in  3  dimensions  when  the  input  data 
are  piecewise  analytic.  This  includes  polyhedron  domains. 

3.  THE  FINITE  ELEMENT  METHOD 


There  are  various  different  forms  of  the  finite  element  method.  We  will 
consider  here  only  the  basic  class  of  finite  element  methods  (for  our  model 
problem) . 

Let  T  -  {r .}  be  a  partition  of  0  into  (in  general  curvilinear)  tri¬ 
angles  or  quadrilaterals  called  elements  r^.  In  the  case  when  T  is  a  tri¬ 
angulation,  we  are  making  the  standard  assumptions.  For  the  general  case  we 
refer  to  [11]  [12],  [36].  We  will  formulate  here  the  assertions  in  the  case 
of  triangulation  only,  although  they  hold  in  general.  Let  H(p,T)  -  [u  € 
H1(fl)|u|T  ,  €  T,  is  a  polynomial  of  degree  p}  be  the  finite  element 

space.  If  is  a  rectangle,  then  polynomials  are  of  degree  p  in  both 

variables.  If  the  elements  are  curvilinear,  then 
"pull-back"  polynomials. 

We  have  assumed  that  the  degree  of  polynomials  are  the  same  over  all  ele¬ 
ments.  The  theory  is  also  developed  for  general  case  when  the  degree  p  can 
be  different  on  different  elements . 


'ri 


are  the  standard 


Let  Hp(p.T)  _  H(p,T)  0  hJ(A)  where_  hJ(0)  -q{u  €  H1(Q)|u  -  0 


on 


r°i 


and  Hq(p,T)  be  the  restriction  of  H(p,T)  on  We  will  assume  that  a 

projection  operator  Pp(p,T)  °f  function  tp  into  Hq(p,T)  be  given  and  we 
denote  o>p  j  "  Pp(p.P)<p* 

The  finite  element  method  consists  now  in  finding  uFg  -  u(p,T)  €  H(p,T) 
so  that 

uFE(p,T)  -<Pp>T  on  r° 


1) 


2) 


U  (air  +  air  fe(‘lxitoi  ■  J,  +*  »•  +  If  *» 


n 


n 


holds  for  any  v  f  Hp(p,T). 

We  will  be  interested  here  in  the  accuracy  of  the  finite  element  solution 
measured  in  the  energy  norm.  Define  e  -  u  -  uPE.  and  let 

,2 


!|e(p  ,T)  I 


ff  ((flf  *  (flf 


FE 

3e 


be  the  error  measured  in  the  energy  norm. 

Two  kind  of  operators  Pp(p,T)  will  be  considered.  Let  y  c  ru  be  the 

side  of  the  triangle  r  with  end >oints  A,B  and  assume  that  y  ■  [-1,1].  A 

-  -1,  B  -  1.  Then  cp_  t!-*  “  i(x)  +  w(x)  where  i(x)  is  linear  function 

P »  «  7 


3 


I 


on  7  such  that  <P  p  l^(+l)  —  <P(-tl)  and  a)  in  the  case  of  H^-projection: 

«  ,T  p-1 

P$<p.T):  <*>'(x)  -  Y  aki/(x),  w'  (x)  -  Y  where  are  the  Legendre  poly- 

k-0  0  1/2 

nomials ;  b)  in  the  case  of  H1'  -projection  Pq  (p,T):  cp(x)  -  Y  a  T,  (x) , 
p  k-0  n  * 

w(x)  =  X  a^T^Cx) ,  where  Tj^x)  are  the  Tchebyschev  polynomials.  Other 
k-0 

projections  are  analyzed  in  [12]. 

As  we  have  seen  there  is  a  large  freedom  in  the  selection  of  H(p,T) 
namely  the  degree  p  and  the  partition  T  and  in  the  selection  of  the 
operators  Pq(P  .T).  We  expect  that  ||e||g  -*  0  if  dim  Hq(p,T)  -*  <®. 

It  is  convenient  to  distinguish  in  this  context  three  versions  of  the 
finite  element  method. 

a)  The  h-version.  Here  a  sequence  (family)  H(p,T^)  is  considered 

when  p  is  fixed  (usually  p  -  1,2)  and  the  mesh  T ^  is  successively  re¬ 
fined  so  that  the  size  h  of  the  elements  of  goes  to  zero. 

b)  The  p-version.  Here  the  mesh  T  is  kept  fixed  and  p  -+  «■>  uniformly 
on  selectively. 

c)  The  h-p  version.  In  this  version  the  mesh  is  simultaneously  refined 
and  the  degree  p  increased  uniformly  or  selectively. 

The  h-version  of  the  finite  element  method  is  the  standard  one.  The 
p-version  is  a  recent  development.  The  first  theoretical  paper  about  the  p- 
version  [28]  and  the  h-p  version  [7]  appeared  in  1981  and  various  results  were 
obtained  since  then. 

There  are  many  codes,  research  and  commercial  utilizing  the  h-version. 

The  only  commercial  code  using  p  and  h-p  versions  is  the  code  PROBE  which  was 
developed  by  NOETIC  Technologies,  St.  Louis,  MO  [49]  [60],  PROBE  solves  two 
dimensional  problems  of  linear  elasticity,  stationary  heat  problems  and 
thermoelasticity  problems  in  two  dimensions  and  elasticity  in  three  dimen¬ 
sions.  Three  dimensional  research  code  STRIPE  was  also  developed  by  Swedish 
Aeronautical  Research  Institute,  see  e.g.  [1]  [2]  [3].  These  codes  have  in 
addition  various  features  as  adaptive  approaches,  various  a-posteriori  error 
estimation,  etc. 

4.  THE  BASIC  PROPERTIES  OF  THE  p  AND  h-p  VERSIONS  OF  THE  FINITE  ELEMENT 
METHOD 


After  1980  the  p  and  h-p  versions  of  the  finite  element  method  was 
studied  in  detail  from  various  point  of  view.  We  will  mention  here  some 
essential  illustrative  results. 

In  one  dimensional  setting  the  versions  were  analyzed  in  detail  in  [34], 
Here,  among  other,  the  optimal  meshes  and  p-distribution  has  been  established 
with  upper  and  lower  bounds  of  the  errors  for  the  three  basic  finite  element 
versions . 

In  two  dimensional  setting  the  following  theorem  is  characteristic  for 
the  performance  of  the  h-p  version.  (For  details  see  [9],  [11],  [12],  [15], 
[35],  [36], [37].) 

Theorem  4.1.  Let  the  solution  u  of  the  problems  (2.1),  (2.2)  belongs 

to  the  set  .  Then  there  is  a  sequence  of  meshes  T ^  and  the  degrees 

p  ^  such  that 


1 

<4 

■ 


» 


4 


t 


(4.1) 


3_ 

-a7N 

II e  II E  <  C  e  ,  a  >  0 


where  -  dim  Hq(p£,T^)  is  the  number  of  degrees -of -freedom  for  the  h-p 
version.  0  □ 


In  one  dimension  the  rate  is 


It  has  been  proven  in  [34]  that 


_  2 

the  optimal  mesh  is  a  geometric  one  with  the  factor  (/ 2  -  1)  =  .17.  The 

experience  shows  that  the  geometric  mesh  with  the  factor  =  .15  is  also 
optimal  in  two  dimension.  Theorem  4.1  together  with  Theorem  2.1  shows  that 
practically  in  any  problem  of  structural  mechanics  the  exponential  rate  of 
convergence  can  be  achieved. 

For  the  p-version  the  following  theorem  is  another  typical  one  (for  more, 
see  [23],  [24],  [25],  [26]). 


Theorem  4.2.  Let  us  consider  the  problem  (2.1)  (2.2)  and  let  /3  - 
max(/J^)  given  in  Theorem  2.1.  Then  for  the  p-version  we  have 

(4.2)  || e || £  < 

while  for  the  h-version  with  uniform  mesh 


(4.3)  1  e  ||  £  >  CN 

In  Theorems  4.1  and  4.2  either  the  projection  Pq  or  Pq^  could  be 
used  provided  <p ,  is  sufficiently  smooth.  The  difference  between  these  two 
operators  occurs  for  the  p-version  when  the  boundary  condition  4)  is  un- 

smooth,  e.g.  tp  €  H®(D,  1/2  <  6  <  3/4.  In  this  case  the  projection 

has  to  be  used.  For  the  h-p  version  there  is  no  difference  in  the  asymptotic 

rate  but  some  difference  occurs  in  the  constant  of  the  estimates.  For  the 
analysis  of  the  influence  of  the  operator  Pq  on  the  accuracy  we  refer  to 
[12]  [15]. 

So  far  we  have  assumed  that  the  domain  0  is  bounded.  Nevertheless  the 
exponential  rate  of  convergence  (4.1)  holds  also  for  the  problem  on  O  - 

-  Q  when  f  has  bounded  support.  Here  the  infinite  elements  and  properly 
selected  shape  functions  have  to  be  used.  For  more,  see  [17]. 

The  h,  p  and  h-p  versions  have  different  aspects  with  respect  to  the  pol¬ 
lution  problem.  In  presence  of  a  singular  behavior  of  the  solution  (e.g.  in 

the  neighborhood  of  the  entrant  corner  of  the  domains)  the  L,,,  error  is  very 

large  in  the  element  consisting  the  comer.  This  effect  disappears  in  ele¬ 
ments  which  are  separated  away  from  the  singularity  by  few  elements.  This 
effect  is  essential  for  a  proper  mesh  design  in  practical  computation.  For 
details  we  refer  to  [16], 

So  far  we  have  dealt  with  the  problem  (2.1)  and  (2.2)  of  second  order. 

For  the  analyses  of  the  finite  element  solution  for  the  problems  of  order  2k 
we  refer  in  the  case  of  the  h-p  version  to  [35]  and  the  p-version  to  [59]. 

For  the  basic  analysis  of  the  p-version  in  3  dimensions,  we  refer  to  [31]  and 
[32].  The  eigenvalue  problem  is,  as  is  well  known,  directly  related  to  the 
"source"  problem  we  addressed  earlier.  See  e.g.  [19]  and  [20],  The  eigen¬ 
functions  belong  to  B^(fi)  and 
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\\-XFE I  <  C  e-^w  ^ 

iUFE(V  '  U<A)!Ie  *  C  ' 

For  more  details,  see  e.g.  [14]. 

5.  THE  IMPLEMENTATION 

There  are  some  essential  features  of  implementation  of  the  p  and  h-p  ver¬ 
sions.  See  e.g.  [60].  The  elements  are  of  a  hierarchical  type  which  leads  to 
augmentation  (bordering)  of  the  local  stiffness  matrices  when  p  is  increased. 
This  also  allows  to  change  very  flexibly  the  degree  of  the  shape  functions 
from  one  element  to  the  other  one.  The  shape  elements  are  (in  two  dimensions) 
of  nodal  type,  side  type  and  internal  type  and  are  based  (in  PROBE)  on  the 
integrals  of  the  Legendre  polynomials .  This  is  important  for  numerical  sta¬ 
bility  aspects. 

The  solvers  can  be  direct  or  iterative.  In  [6]  we  have  discussed  the 
complexity  of  these  solvers  in  the  framework  of  parallel  computation  in  the 
dependence  on  the  architecture  of  the  computer.  In  [48]  we  have  proven  vari¬ 
ous  aspects  of  the  preconditional  conjugate  gradient  method.  It  has  been 
shown  that  is  is  essential  to  eliminate  first  all  element- internal  unknowns 
—  which  can  be  done  completely  in  parallel  —  and  then  for  properly  designed 

2 

conjugate  gradient  method  the  number  of  iterations  depend  only  on  lg  p  (p 
degrees  of  elements). 

6.  THE  PROBLEM  OF  THE  MESH  DESIGN 

One  of  the  most  laborious  part  of  the  finite  element  analysis,  especially 
in  three  dimensions,  is  the  mesh  generation  also  when  sophisticated  mesh 
generators  are  used.  The  use  of  large  elements  (possibly  of  high  degrees) 
which  are  describing  only  the  geometry,  greatly  simplifies  the  user’s  work, 
also  if  possibly  on  expense  of  the  computer  time.  (It  is  necessary  to  realize 
that  the  relation  between  manpower  cost  and  computation  cost  is  going  steadily 
up.)  The  option  of  a  change  of  the  degrees  of  elements  increases  signifi¬ 
cantly  the  flexibility  of  the  program  and  gives  the  user  effective  tool  for 
the  quality  control.  The  p  and  h-p  version  programs  give  such  options.  It  is 
advantageous  therefore  to  create  directly  or  indirectly  the  proper  mesh  and  to 
achieve  then  the  desired  accuracy  by  an  appropriate  choice  of  the  degrees, 
which  can  be  made  in  an  adaptive  mode.  The  goal  is  to  achieve  the  same  combi¬ 
nation  of  the  degrees  and  mesh  refinement  which  would  be  obtained  for  the 
given  accuracy  by  the  h-p  version  directly.  To  achieve  this  goal  two  avenues 
could  be  followed,  the  expert  system  and  the  adaptive  approach.  The  expert 
system,  see  e.g.  [21],  [51],  advises  the  user  how  to  design  the  mesh  and 
element  degrees  for  the  requested  accuracy  and  provides  the  user  with  a  mesh 
generator . 

The  adaptive  approach  (see  e.g.  [1],  [2],  [3],  [34]  Part  3)  which  is  pos¬ 
sible  to  see  as  an  "automatic"  expert  system  makes  various  decisions  for  the 
users.  Both  approaches  have  some  common  parts  but  the  concepts  are  signifi¬ 
cantly  different.  We  refer  also  to  [46],  [52],  [56],  [65],  [66],  [67]  for 
various  additional  aspects . 
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7.  THE  ROBUSTNESS 


An  effective  method  has  to  perform  uniformly  well  for  a  broad  class  of 
input  data.  The  elasticity  problems  can  be  in  practice  nearly  singular  as, 
for  example,  in  the  case  of  nearly  incompressible  material,  various  plate  and 
shell  theories,  in  the  case  of  thin  domains,  etc.  The  h-version  suffers  in 
these  cases  by  the  "locking"  problem  which  has  to  be  overcome  by  various 
special  approaches  as  reduced  integration,  etc.  Problems  of  these  types  are 
avoided  when  the  p  and  h-p  versions  are  used.  The  convergence  rate  in  the 
energy  norm  then  (in  contrast  to  the  h-versions)  is  uniform  with  respect  to 
the  Poisson  ratio  when  higher  degrees  of  elements  are  used.  See  e.g.  [56], 
[39],  [73].  For  computation  of  the  pressure  for  v  ~  1/2,  see  [69], 

8.  THE  QUALITY  CONTROL  OF  THE  SOLUTION 

It  is  essential  to  have  a  possibility  of  a  quantitative  assessment  of  the 
quality  of  computed  data.  For  the  survey  of  today's  general  ideas  and  results 
in  this  direction,  we  refer  to  [45].  In  the  case  of  the  p  and  h-p  versions 
there  is  relatively  easy  way  for  the  quality  assessment  of  any  data  of  the 
interest  by  changing  the  degrees  and  by  an  extrapolation  procedure.  This  ap¬ 
proach  is  very  effective  because  it  indicates  reliably  the  errors  of  any  com¬ 
puted  data  of  interest,  the  energy  norm,  value  of  stresses,  stress  intensity 
factor,  etc.  See  also  [1],  [11], [33],  [36],  [38],  [49],  [58],  [60],  [61], 
[64],  [68]. 

9.  p  AND  h-p  VERSION  IN  OTHER  PROBLEMS 

So  far  we  mentioned  only  the  use  in  the  context  of  elliptic  problems 
[57],  [77]  addresses  the  use  of  p-version  for  the  boundary  element  method, 

[74],  [75],  [76]  in  the  context  of  the  Kirchhoff  plate  analysis.  [53],  [71], 
show  that  the  p  and  h-p  versions  are  very  effective  tools  for  computation  of 
thin  structures  as  plates  and  shells  and  [54]  for  optimal  design.  [13] 
analyzes  the  p-version  for  parabolic  equations  when  high  order  polynomials 
are  used  in  space  and  time  and  [41]  is  using  the  p-version  for  solving 
elliptic  stochastic  equations. 

10.  THE  COMPUTATIONAL  AND  ENGINEERING  EXPERIENCE 


Because  of  the  developed  commercial  code  PROBE  and  research  code  STRIPE 
an  extensive  experience  is  already  available  in  the  research  and  industrial  < 

use.  For  the  industrial  experience  we  refer  e.g.  to  [29],  [49]  where  numer¬ 
ical  results  are  presented.  For  the  numerical  data  we  refer  e.g.  to  [1], 

[2],  [3],  [4],  [11],  [12],  [13],  [14],  [15],  [16],  [17],  [21],  [22],  [36], 

[38],  [41],  [48],  [51],  [53],  [54]  [61],  [62], [63],  [66],  [67],  [68],  [71], 

[72].  The  experience  shows  that  the  p  and  h-p  versions  of  the  finite  element 

method  has  many  practical  advantages  for  the  engineering  computations  for  < 

linear  elliptic  problems. 

11.  RELATION  TO  SOME  OTHER  METHODS 


The  ideas  of  the  p  and  h-p  versions  are  related  for  example  to  various 
methods  used  in  fluid  dynamics  as  the  spectral  method  and  its  variation.  We 
refer  e.g.  to  [46]  and  references  given  there.  A  commercial  fluid  dynamics 
code  NEKTON  which  is  based  on  spectral  method  and  developed  by  NEKTONICS  INC, 
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Badford,  MA  is  closely  related  to  the  idea  of  h-p  version  of  the  finite 

element  method.  See  also  [30]. 
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